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Abstract: We present phenomenological results for vector boson pair production at the 
LHC, obtained using the parton- level next-to- leading order program MCFM. We include 
the implementation of a new process in the code, pp — t- 77, and important updates to 
existing processes. We incorporate fragmentation contributions in order to allow for the 
experimental isolation of photons in 77, W^, and Z7 production and also account for 
gluon-gluon initial state contributions for all relevant processes. We present results for a 
variety of phenomenological scenarios, at the current operating energy of y/s = 7 TeV and 
for the ultimate machine goal, = 14 TeV. We investigate the impact of our predictions 
on several important distributions that enter into searches for new physics at the LHC. 
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Figure 1: NLO boson production in pp-collisions. The decay branching ratios of the VF's and Z's 
into one species of leptons are included. For 77 and V'f we apply pr cuts of 25 and 10 GeV to 
photons respectively. 

1. Introduction 

The current plan for the LHC calls for running in both 2011 and 2012. Running in 2011 
is at a centre of mass energy at ^/s = 7 TeV, with a baseline expectation of 1 fb~^ per 
experiment and a good chance that greater luminosity will be accumulated. At the end of 
the 2012 run it is likely that data samples in excess of 5fb~^ will have been accumulated 
by both of the general purpose detectors. Data samples of this size will (at the very least) 
allow detailed studies of the production of pairs of vector bosons. 

It therefore seems opportune to provide up-to-date predictions for the production of 
all pairs of vector bosons, specifically for the LHC operating at 7 TeV. This extends the 
previous implementation of diboson production in MCFM [1] which was focussed primarily 
on the Tevatron. Moreover, we also consider the production of final states that contain real 
photons. This requires the inclusion of fragmentation contributions in order to address the 
issue of isolation in an experimental context. In addition, we have also included the con- 
tribution of the gluon-gluon initial state to a number of processes. These finite corrections 
are formally of higher order but can be of phenomenological relevance at the LHC where 
the gluon flux is substantial. 

A review of the current experimental status of vector pair boson production, primarily 
from the Tevatron, can be found in ref. [2]. The production of pairs of vector bosons is 
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crucial both in order to check the gauge structure of the Standard Model (SM) and in the 
search for new physics. This is because production of vector boson pairs and the associated 
particles from their decay, enter as irreducible backgrounds for many Higgs and new physics 
searches. The observationally most promising decays of the Higgs boson are to two photons 
(for a light Higgs), or to two W^s or two Z's that decay leptonically. Clearly vector boson 
pair production is an irreducible background in these searches. Processes with leptons and 
missing energy are typical signatures of many new physics models, of which supersymmetry 
is a classic example. Again, knowledge of SM processes which possess multiple leptons and 
missing energy is crucial in the quest to discover or rule out these models. 

In Fig. |l] we show the rates for various electroweak processes at energies between y/s = 7 
and 14 TeV. This figure serves both as a road-map to this paper and as an indication of 
the relative size of the various diboson processes. We present the cross sections for single 
boson production to illustrate the orders of magnitude which separate single boson and 
diboson production. Where appropriate we have included the branching ratios of vector 
bosons to a single family of leptons and applied a transverse momentum cut of 10 GeV 
{W^ and Z7) and 25 GeV (77) to photons. No other cuts are applied to the boson decay 
products. 

Updating the diboson processes in MCFM for the new energy range probed at the 
LHC is the primary aim of this work. With that in mind we begin in section |2| by outlining 
the steps needed to include photon fragmentation in the code. Section ^ serves as an 
overview, describing the parameters that we use and outlining the processes that receive 
extra corrections from gluon initiated production mechanisms. Section ^ discusses the 
phenomenology of 77 production at the LHC. We investigate the role of isolation on the 
cross section and the impact of Higgs search cuts on di-photon production. Sections ^ and ^ 
contain our predictions for 1^7 and Z7 production at the LHC. We investigate the role 
of final-state radiation in our calculations and compare our NLO results with the recently 
reported cross sections from CMS [3]. Sections 0, ^ and ^ turn to the production of two 
massive vector bosons. We are able to compare our prediction for the WW cross section 
with early results from ATLAS and CMS [4, 5]. We examine the effect of the gluon initiated 
processes in the WW and ZZ final states, with particular emphasis on their role as Higgs 
backgrounds. For WZ production we discuss briefly the properties of boosted Z's. Finally 
in section |l^ we draw our conclusions. Appendix |A| contains a more detailed discussion 
of our electroweak parameters whilst appendix ^ presents formulae for the gg — )• V1V2 
amplitudes as implemented in MCFM. 

2. Photon fragmentation 

Since we will consider a number of final states including photons we must first discuss 
the additional complications that this involves, compared to the production of W and Z 
bosons. Experimentally, the production of photons occurs via two mechanisms. Prompt 
photons are produced in hard scattering processes whilst secondary photons arise from the 
decays of particles such as the vr'^. Since secondary photons are typically associated with 
hadronic activity one can attempt to separate these contributions by limiting the amount 
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of hadronic energy in a cone of size Rq = ■y/(A?7^ + A(p'^) around the photon. Experimental 
isolation cuts are of the form, 



^T(had) < ehpl or ^Erihad) < E^""" . (2.1) 



eRo eRo 



Thus the transverse hadronic energy, Et (had) , is limited to be some small fraction of the 
transverse momentum of the photon or cut off at a fixed, small upper limit. 

Matters are complicated both experimentally and theoretically by a second source of 
prompt photons. A hard QCD parton can fragment non-perturbatively into a photon. As 
a result a typical photon production cross section takes the form, 



a^{Ml) + j dzD''{z)a''{z,Ml). (2.2) 



Here a'^ represents the direct component of the photon production cross section whilst the 
second term arises from the fragmentation of a parton a into a photon with momentum zp"". 
Each contribution separately depends on the fragmentation scale, Mp- The fragmentation 
functions, taken as solutions to a DGLAP equation are of (leading) order uew /o-s- This 
means that they are formally of the same order as the leading order direct term. At high- 
energy hadron colliders, the QCD tree-level matrix element, coupled to a fragmentation 
function can become the dominant source of prompt photon production. However, the 
magnitude of these terms can be drastically reduced by applying the isolation cuts described 
above. This is due to the fact that the fragmentation functions strongly favour the low 
z region. Once the photon is isolated, z is typically large enough that the fragmentation 
contribution drops substantially from the unisolated case. 

A theoretical description of isolated photons is complicated because of the occurence 
of collinear singularities between photons and final-state quarks. A finite cross section 
is only obtained when these singularities are absorbed into the fragmentation functions. 
As a result the only theoretically well-defined NLO quantity is the sum of the direct and 
fragmentation contributions. Once these two contributions are included one can isolate the 



photon using the cuts of Eq. (2T) in an infrared safe way [6]. 

Although the underlying dynamics of photon fragmentation are non-perturbative the 
evolution of the functions with the scale Mp is perturbative. In the same manner as 
the parton distribution functions, the fragmentation functions satisfy a DGLAP evolution 
equation. In MCFM we use the fragmentation functions of ref. [7], which are NLL solutions 
to the DGLAP equation. 

Final state quark-photon collinear singularities are removed using a variant [6] of the 
Catani-Seymour dipole subtraction formalism [8]. More specifically, we treat the photon in 
the same manner as one would treat an identified final state parton (with the appropriate 
change of colour and coupling factors). Integration of these subtraction terms over the 
additional parton phase space yields pole pieces of the form [6], 



1 r(l-e) / 4V \ a 2 
e r(l - 2e) V Ml ) 27r ''''' 
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where P^q{z) is the tree level photon-quark splitting function. This piece is the lowest 
order definition of the photon fragmentation function in the MS scheme. This singularity 
is then absorbed into the fragmentation functions to yield finite cross sections. 

Since the isolation cuts reduce the magnitude of the fragmentation contributions we 
calculate the QCD matrix elements a"'{z,Mp) to LO, i.e. we neglect NLO corrections to 
the fragmentation processes. 

An alternative procedure, in which one can avoid calculating the fragmentation con- 
tributions altogether, is to follow the smooth cone isolation of Frixione [9]. In such an 
approach one applies the following isolation prescription to the photon, 

Using this prescription, soft radiation is allowed inside the photon cone but collinear sin- 
gularities are removed. Since the smooth-cone isolation is infra-red finite, there is no need 
to include fragmentation contributions in this case. Currently this isolation is difficult to 
implement experimentally and therefore it is not used in this paper. ^ 



3. Overview 

The results presented in this paper are obtained with the latest version of the MCFM code 
(v6.0). We use the default set of electroweak parameters as described in Appendix]^. 

For the parton distribution functions (pdfs) we use the sets of Martin, Stirling, Thorne 
and Watt [10]. For the calculation of the LO results presented here we employ the corre- 
sponding LO pdf fit, with 1-loop running of the strong coupling and as{Mz) = 0.13939. 
Similarly, at NLO we use the NLO pdf fit, with as{Mz) = 0.12018 and 2-loop running. 
The fragmentation of partons into photons uses the parametrization "set 11" of Bourhis, 
Fontannaz and Guillet [7]. 

As mentioned in the introduction, for several processes we have included contributions 
of the form gg — )• V1V2. These contributions proceed through a closed fermion loop and 
form a gauge invariant subset of the one-loop amplitudes. However, since there is no gg 
tree level contribution the first time these pieces enter in the perturbative expansion is at 
a|, (i.e. NNLO). Simple power counting would thus lead one to assume that these pieces 
are small, of the order of a few percent of the LO cross section. At the LHC this is often 
not the case, since the large gluon flux in the pdfs can overcome the O^Og) suppression in 
the perturbative expansion. The resulting gluon-gluon contributions are instead 0(10%) 
of the LO cross section, i.e. these pieces are comparable to the other NLO contributions. 

Charge conservation ensures that not all diboson processes receive these gluon-gluon 
initiated contributions. The allowed processes are gg — )• {'yy, Z^,W~^W~ , ZZ}, each of 
which has been studied in some format in the past [11-26]. We refer the reader to the appro- 
priate section for details of each calculation. We note that for gg — )• {Zj, W^W^ , ZZ} — )• 

^Smooth cone isolation is however available in MCFM for theoretical comparisons. 



- 5 - 



leptons, we present (to the best of our knowledge) analytic formulae for the helicity am- 
plitudes for the first time. These formulae were readily obtained from the amplitudes for 
the process e"'"e~ — >• 4 partons [27]. 

Since there are no gg tree level contributions, each of these 1-loop amplitudes is both 
infrared and ultraviolet finite. This means that, once calculated, these contributions are 
simple to implement in MCFM. Throughout this paper we include these pieces in the NLO 
results except for the gg — )• 77 section. Since the strong corrections to the gg — )• 77 process, 
including the two- loop amplitude [13] are known and sizeable [14], we proceed differently 
for this process. The two-loop gg — )• 77 amplitude is infrared divergent and must be 
combined with real radiative corrections in exactly the same manner as a canonical NLO 
contribution [14]. As a result these contributions are included in our NLO predictions for 
the diphoton process, while the one-loop gg — )• 77 calculation is included in the LO result. 

4. 77 production 

4.1 Description of the calculation 

In view of its role as the principal background in the search for the light Higgs boson 
in the decay mode -ff — t- 77, it is important that the prediction for Standard Model 
diphoton production is as accurate as possible. The production of photons in hadron- 
hadron interactions proceeds through the Born level process, 

q + (4.1) 

Corrections to this picture due to QCD interactions have been first considered at 0{as) in 
ref. [28] and the results for that process have been included in the Diphox Monte Carlo [29]. 
The large flux of gluons at high energy - in particular at current LHC energies - means that 
diagrams involving loops of quarks can give a significant additional contribution [11, 12, 15], 

5 + 5^77- (4.2) 

Since these contributions can be rather large, in order to obtain a reliable estimate of 
their contribution to the diphoton cross section it is necessary to include higher order 
corrections. The results of such a calculation, involving two-loop virtual contributions [13], 
were presented in ref. [14]. 

The results presented in this section are obtained using our current implementation 
in MCFM which is as follows. The gg process is included at NLO using the two loop 
matrix elements of ref. [13] and following the implementation of ref. [14]. We include 
five flavours of massless quarks and neglect the effect of the top quark loops, which are 
suppressed by 1/mf. Next-to-leading order corrections to the qq initiated process are more 
straightforward to include, although some care is required due to the issues of photon 
fragmentation and isolation that have been described in section 

We can compare our implementation ofpp — t- 77 to Diphox [29]. Diphox contains NLO 
predictions for both the direct and fragmentation pieces, but only includes the gg initiated 
pieces at leading order. In MCFM we include NLO predictions for the direct pieces, LO 
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[TeV] 



7 



10 
11 
12 
13 
14 



a^O(77) [Pb] 



35.98(0) 
43.04(1) 
50.32(1) 
57.76(1) 
65.37(1) 
73.07(1) 
80.89(1) 
88.76(2) 



NLO 



(77) [pb] 



1+5% 

'-6% 

1+4% 
'-6% 

+5% 
-5% 



73.0(2)+^^° 



47.0(1 



55.8(1 



64.3(1) 



81.8(2) 



5% 

+3% 
-5% 



90.5(3)+^J 

99.1(3)+^^" 

1-3% 

-5% 



108.1(3)+^^" 



Table 1: L O and NLO cross sections for diphoton production at the LHC with the acceptance cuts 
of Eq. (4.3), as a function of i/s. The Monte Carfo integration error on each prediction is shown 
in parentheses. For the NLO results the theoretical scale uncertainty is computed according to the 
procedure described in the text and is shown as a percentage deviation. 



predictions for the fragmentation processes (using NLL fragmentation functions) and the 
"NLO" gg predictions. For isolated photons the "NLO" gg corrections represent around 
5% of the total cross section, so we expect them to be at least as important as the NLO 
corrections to the fragmentation piece. 

4.2 Results 

As a point of reference, we first consider the cross section for unisolated photons at the 
LHC, for various centre-of-mass energies. We apply only basic acceptance cuts on the two 
photons, 

plj. > 25 GeV , |r/^| < 5 . (4.3) 

The cross sections we report are completely inclusive in any additional parton radiation. 
For our theoretical predictions we choose renormalisation (/x/j), factorisation {^jlf) and 
fragmentation scales {Mp) all equal to the diphoton invariant mass, m^^. The results of 
our study at LO and NLO are shown in Table ||, where the percentage uncertainties quoted 
on the NLO cross sections are estimated by varying all scales simultaneously by a factor of 
two in each direction. The inclusion of both gg and qq processes in the LO result, and 
the next order corrections to both at NLO, results in only a mild 20-30% increase in the 
cross section at NLO. Moreover, these predictions are rather stable with respect to scale 
variations over the range studied, with deviations in each direction of at most 6%. 

We now wish to investigate a more realistic set of cuts in which the photon is iso- 
lated. Since this final state is particularly interesting in the context of a low-mass Higgs 
search [25], for illustration we adopt the set of cuts used in an early search by the AT- 
LAS collaboration [30]. The photons are required to be relatively central and subject to 
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Figure 2: The NLO prediction for the diphoton cross section (in picobarns) as a function of the 
cent re of mass energy, -^s. The cross sections are shown for three sets of cuts: only the basic cuts of 
Eq. (O) (upper, blue curve); the staggered cuts of Eq. (14) (middle, magenta curve); the isolated 
photon cross section, Eqs. (4.4, E^) (lower, red curve). 



staggered transverse momentum cuts, 

p^i>40GeV, p:^'>25GeV, |r/^J < 2.5 , (4.4) 

and are isolated using a fixed maximum hadronic energy in a photon cone (c.f. Eq. ( |2.1| )), 

Ro = 0.4 , = 3 GeV . (4.5) 

The effect of these cuts, as a function of y/s, is shown in Fig. The effect of the staggered 
cuts, Eq. (|4.4| ), is to lower the cross section by approximately a factor of three compared 
to the basic cuts of Eq. (|4.3|) . The isolation condition, Eq. ( |4.5| ), further reduces the cross 
section from the nominal unisolated prediction by about 9%. We note that this reduction 
is smaller than one would typically expect when going from unisolated to isolated cross 
sections. This is due mostly to the staggered cuts which favour the 3 particle final state. 

In fact the cross section is rather insensitive to the amount of transverse hadronic 
energy allowed in the isolation cone. This is illustrated in Fig. ^ which shows the de- 
pendence of the cross section on the value of the isolation parameter E^^^. As a result 
of the small variation over this range, isolation cuts of the form E + dp^p where E and 5 
are constants and 6 <^ E are well-approximated theoretically by using a simple constant 
= E + 5pl^^. 

For the cross sections presented so far we have chosen to set all scales entering our 
calculation equal to the invariant mass of the two photons, = fip = Mp = /iO) with 
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Figure 3: The fraction of tlic unisolated diplioton cross section tliat remains when the photon is 
isolated, as a function of the maximum amount of transverse hadronic energy allowed in the photon 
isolation cone, E^^^. The centre of mass energy is y/s = 7 TeV and photons are identified according 
to the staggered cuts of Eq. (4.4). The radius of the isolation cone is Rq = 0.4. 



fiQ = m^^. In order to illustrate the impact of this choice, in Fig. Q we show the dependence 
of the theoretical predictions on the common scale fi when it is varied by a factor of four 
about ^0- In addition to the scale dependence of the total predictions we also consider the 
variation with ^ of the individual partonic channels that appear at each order. Although the 
scale dependence of the individual partonic processes is typically quite large (for example, 
for the qq and qg initiated processes at NLO), the sum over all contributions is relatively 
scale-independent. The LO cross section in particular has a tiny variation in this range. 
The fact that the NLO corrections are large and not reproducible for any choice of scale 
considered at LO serves as a reminder that the scale variation is not indicative of the 
theoretical uncertainty at that order. We also note the large K-factor when going from LO 
to NLO (~ 3.2), which is in stark contrast to the mild corrections observed when imposing 
only basic acceptance cuts (c.f. Table ||). This difference can easily be understood from 
the nature of the cuts in Eqs. ( [4.4| , [1.5| ). For the Born and virtual contributions the 
staggered pT cut is effectively a pj? > 40 GeV cut due to the 2 — )• 2 kinematics. Photons of 
Pt < 40 GeV can only be produced by fragmentation or real radiation diagrams in which 
a parton is available to balance the staggered transverse momenta. As a result these cuts 
strongly favour real radiation diagrams, a fact that is also evident from the size of the qg 
contribution in Fig. g Thus we would advocate using equal cuts on the photons, since 
they do not emphasize the role of higher order corrections. 



Finally, we consider predictions for the diphoton invariant mass distribution, a key 
ingredient in the search for a light Higgs boson. Our results for y/s = 7 TeV and the cuts 
of Eqs. {iA, are shown in Fig. It is clear that, in order to provide a good prediction 
for this distribution, one must include not only the gluon-gluon initiated process but also 
the NLO corrections throughout. 
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Figure 4: Dependence of the LO and NLO diphoton cross sections at y/s — 7 TeV (in pb) on the 
scale choice fj,. We vary fi = ii^j — ~ Mp about the central scale choice /io — rn^-y. Total cross 
sections are shown in black whilst colours are used to denote the scale dependence of particular initial 
states: quark-antiquark (red), quark-gluon (magenta), gluon-gluon (blue). Photons are defined and 
isolated according to Eqs. (4.4, 4.5). 



5. W 'y production 

5.1 Description of the calculation 

The production of a W boson and a photon proceeds at Born level via quark-antiquark 
annihilation, 

q + q'^W'^'y. (5.1) 

This process was first calculated several decades ago [31], with the effect of radiative correc- 
tions subsequently accounted for in ref. [32]. Since then the subject has been revisited sev- 
eral times. A fully differential Monte Carlo implementation of the NLO result is presented 
in ref. [33] , making use of the helicity amplitudes calculated in ref. [34] . Spin correlations 
in the decay of the W boson are included although no photon radiation from the lepton 
is allowed. Electroweak corrections to this process [35] and NLO QCD corrections to the 
related W^7+jet final state have also been computed [36]. 

In this section we present results using the current implementation of this process in 
MCFM. The diagrams that contribute to this process at leading order are shown in Fig. |6|. 
The next-to-leading order diagrams are obtained by dressing these diagrams with both 
virtual and real gluon radiation. The contribution to the full amplitude arising from three 
of these diagrams is readily obtained from the helicity amplitudes of ref. [34]. The final 
diagram, including appropriate dressings that are straightforward to compute, accounts for 
the additional contribution from photon radiation in the leptonic decay of the W boson. 
The resulting amplitude retains full spin correlations in the decay. 
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Figure 5: The diphoton invariant mass distribution at -^s = 7 TeV (in fb/GeV). We apply the 
staggered cuts described in the text and indicate LO results with dashed curves and NLO results 
with solid curves. The two upper (blue) curves show the full predictions at a given order, while the 
lower (red) curves indicate the gluon-gluon initiated contributions only. 




(a) 





(c) (d) 

Figure 6: Leading order diagrams for Iv)^ production. The diagrams (a),(b) and (c) can 

be considered as radiation in the production process, while the final diagram (d) corresponds to 
photon radiation from the lepton in the W decay. 



5.2 Results 

In order to define the final state for this process we apply a basic set of kinematic cuts, 

V^j, > 10 GeV , Re^ > 0.7 , (5.2) 



[TeV] 


a^O{e+u-f) [pb] 


^JVLO(g+^^) [pb] 


£7^0(e-I77) [pb] 


^iVLO(g-^^) [pb] 


7 
8 
9 
10 
11 
12 
13 
14 


23.02(6) 
26.86(8) 
30.62(8) 
34.6(1) 
38.4(1) 
42.2(1) 
45.9(1) 
49.8(1) 


30.1(l)+^5 
35.1(2)+?% 

39.6(2)+5,5 
44.2(4)+^% 

48.8(3)lg 
54.0(4)+3% 
57.7(4)+3% 
62.8(4)+^% 


15.46(5) 
18.53(7) 
21.26(8) 
24.13(8) 
27.1(1) 
30.2(1) 
33.1(1) 
36.0(1) 


21-1(1)^S 
24.6(l)+g° 

28.4(2)+^5 
32.2(2)+3% 

35.7(2)1^5 
39.4(2)1^5 
43.6(3)lg: 
47.4(3)lg 



Table 2: Cross sections for W{-^ £v)j production as a function of energy, using only the cuts of 
Eq. ( ^.5[ ). The cross sections are calculated including the effects of photon radiation in the W decay 
and the central values are obtained using /i/j — fip — Alp = Mw. The uncertainty is derived from 
the scale dependence, as described in the text. 



and demand that the photon be isolated as before, i?o = 0.4 and E'^^^ = 3. In this 
subsection we consider W bosons which decay leptonically. We do not apply any cuts to 
the leptons, except for the photon-lepton separation cut which ensures that the photon- 
lepton collinear singularity is avoided. The resulting cross sections are given, as a function 
of -y/s, in Table |2[ We present results for the LO and NLO cross sections for €"''1^7 and 
e~T''y separately. The cross sections have been calculated using a central scale choice of 
= fJ'F = Mp = M\Y, with upper and lower extrema obtained by evaluating the cross 
section at {hr = M]v/2,fj,p = 2M\y} and {fiR = 2Miv,fip = Mw/2} respectively. The 
fragmentation scale is kept fixed at M\y throughout since its variation does not lead to a 
significant change in our results over the range of interest. Prom this table we can readily 
extract our NLO prediction for the W'j cross section (summed over both W+ and W) at 
current LHC operating energies with the cuts and isolation described above, 

a^^^ipp ^Wj + X)x BR{W eu) = 51.2 tH pb . (5.3) 

This is to be compared with a recently-reported cross section from the CMS collabora- 
tion [3]. They find, 

a'^'^^^ {pp ^ Wj + X) X BR{W (.v) = 55.9±5.0 (stat)±5.0 (sys)±6.1 (lumi) pb , (5.4) 

in good agreement with the Standard Model expectation. 

Varying the factorisation and renormalisation scales in the manner that we have chosen 
requires some further justification since the normal theoretical preference is to change them 
together in the same direction. A comparison of these two choices, varying ^p and /x/j in 
the same and opposite directions, is shown in Fig. |^. We observe that there is essentially 
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Figure 7: Scale variation for W~^{-^ e^v)^ production, applying only the basic cuts of Eq. (5.5) 



For the red curve we vary the factorisation and renormalisation functions in opposite directions, 
whilst for the blue curve we vary them in the same direction. The fragmentation scale is kept fixed 
at Mw 



no change in the NLO e^z^7 cross section when the scales are varied in the same direction, 
which is due to the qualitatively different behaviour of the contributing partonic states. 
The qq initial state is dominated by variations in the factorisation scale and grows with 
increasing iip- Conversely, the gq initial state depends most strongly on the renormalisation 
scale and decreases with increasing jjLR. The combination of these two initial states results 
in a very small net scale dependence. Since this is simply a fortuitous cancellation and 
higher order corrections to the NLO cross section will likely not be bracketed by this small 
scale variation, we choose to vary the scales in opposite directions instead. We believe that 
this results in a more credible estimate of the theoretical uncertainty of the calculation. 

At LHC centre of mass energies the dominant contributions to the cross sections that we 
have presented so far result from the radiation of a photon from the lepton in the W decay. 
For studies of anomalous couplings of vector bosons to photons, and for the observation 
of radiation zeros in rapidity distributions, it is most useful to suppress this contribution. 
This is achieved by applying a cut on the transverse mass [Mt] of the photon-lepton-MET 
system, Mt > 90 GeV. To investigate the role of lepton cuts on the cross section and 
distributions we will first present results for the cross section at ^/s = 7 TeV including the 
W decay for various sets of cuts. The three sets of cuts that we will consider are, 

Basic Photon : > 10 GeV, \r]^\ < 5, Ri^ > 0.7, Rq = 0.4, E^"'' = 3 GeV. (5.5) 
Mt cut : Basic Photon + Mt > 90 GeV. (5.6) 
Lepton cuts : Mt cut + > 25 GeV, > 20 GeV, \r]e\ < 2.5. (5.7) 

For each of these sets of cuts we will perform our NLO calculation in two different ways. 
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Decay 


Cuts 




cT^^O(e+z/7) 






No FSR 


Basic 7 


4.88 


8.74 


3.15 


6.01 




Mt cut 


1.99 


3.78 


1.26 


2.66 




Lepton cuts 


1.49 


2.73 


0.86 


1.77 


Full 


Basic 7 


23.0 


30.1 


15.5 


21.1 




Mt cut 


2.12 


3.94 


1.34 


2.75 




Lepton cuts 


1.58 


2.85 


0.91 


1.81 



Table 3: W{-^ Iv)^ cross sections in picobarns at ^/s — 7 TeV for the various scenarios described 
in detail in the text. Results in the upper half ("No FSR") correspond to neglecting diagrams 
containing photon radiation in the W decay, while the cross sectio ns i n the lower half ("Full") 
include this effect. The cuts on the final state are specified in Eqs. (|5.5| )-(5.7). Statistical errors 
are ±1 in the final digit. 



In the first case (no final-state radiation, "No FSR") we will omit diagram Fig. |6|(d) corre- 
sponding to photon radiation in the W decay (and its appropriate NLO dressings). Such 
an approach is natural if one demands that the lepton-neutrino system is produced exactly 
on the W mass-shell. This corresponds to the approach taken in ref. [33]. This constraint 
cannot be implemented physically. For the second case ("Full") we follow our usual pro- 
cedure and include this diagram and NLO counterparts. The results are summarised in 
Table |. 

We observe that with just the basic cuts of Eq. (5.5) the difference in predicted cross 
section between the two calculations is very large. At NLO the full result is over three times 
larger than the "No FSR" equivalent. As claimed earlier, applying the Mt cut of Eq. ( ^.6| ) 
significantly reduces this difference. The NLO cross section including radiation in the decay 
is about 3% higher. The quantity that is most relevant experimentally corresponds to the 
full cuts given in Eq. (5/7). In that case the two calculations differ by at most 4% at NLO. 

Including the final state radiation of photons not only significantly increases W'j cross 
sections at the LHC, it also changes the character of the radiation zero [37] that is present 
in the amplitude. The signature of the radiation amplitude zero can be seen in the distri- 
bution of the pseudorapidity difference between the charged lepton and the photon. Our 
predictions for this distribution, using fiji = fip = Mp = M]y and applying the full lepton 
cuts of Eq. (5.7), are shown in Fig. ^. The dashed blue curve in Fig. ^ represents the NLO 
rapidity difference with lepton cuts (Eq. ( |5.7D ), but with no cut on Mt applied. We observe 
that the characteristic dip associated with the radiation zero has been completely filled in 
by the radiation of photons from the charged lepton. This is due to the fact that this 
configuration favours a collinear electron-photon pair so the rapidity difference between 
the two is usually small. Applying the Mt cut (black curve) removes the majority of these 
configurations and the dip is restored. With the Mt cut the NLO prediction from the full 
theory is similar to the result from the "No FSR" calculation (red curve). 
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Figure 8: NLO Predictions for the pseudorapidity difference between the charged lepton and the 
photon in W{—>- iv )^ e vents, for three different levels of the calculation. For all curves we apply the 



lepton cuts of Eq. (5/7). The black curve represents the complete NLO prediction. The red dashed 
curve represents the NLO prediction in the case where no photon radiation is allowed from the 
lepton ("No FSR"). The blue dashed curve has no cut on Mt, but keeps the cuts on the leptons. 



6. Z7 production 

6.1 Description of the calculation 

The production of a Z boson and a photon primarily occurs through the Born process, 

q + (6.1) 

The next-to- leading order corrections to this were computed in refs. [32, 38] and later 
extended to the case of a decaying Z boson in ref. [33]. Electroweak corrections to this 
process have also been computed [35, 39]. 

A further contribution arises from the process, 

g + g^Z^, (6.2) 

which proceeds via a quark loop. Since this contribution is finite it can be computed 
separately, as first detailed in refs. [15, 16]. More recently this process has been computed 
including the leptonic decay of the Z boson and other higher order contributions [17]. 

The results presented in this section are obtained using our current implementation 
in MCFM which is as follows. Strong corrections to the qq initiated process are fully 
included, also allowing additional contributions from fragmentation processes. The gg 
process is included for five flavours of massless quarks. The contribution from massive top 
quark loops is suppressed by 1/mf and is therefore neglected. We find agreement with the 
large top-mass limit of the results presented in ref. [16], where the full dependence on the 
top and bottom quark masses has been kept. 
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Analytic expressions for the gg — >• Z7 amphtudes that we include may be obtained 
from existing results for e~^e~ — )• 4 partons [27], as described in Appendix ^ 

6.2 Results 

We begin by assessing the impact of radiation in the decay of the Z boson to charged 
leptons. As before we consider three sets of cuts to illustrate the difference between the 
two calculations. These are: 

Basic Photon : 171^+^- > 50 GeV, > 10 GeV, |r/^| < 5, Re^ > 0.7, 

Ro = 0.4, = 3 GeV . (6.3) 

Mu^ cut : Basic Photon + > 100 GeV . (6.4) 

Lepton cuts : cut + > 20 GeV, \r]i\ < 2.5 . (6.5) 

The first set of cuts, Eq. (|6^), is very similar to the basic cuts for the process (Eq. ( |5.5| )) 
but with an additional dilepton invariant mass cut in order to select real Z events. The 
cut on the transverse mass Mp has been replaced with a cut on the invariant mass of 
the photon+leptons system. This reflects the fact that for Z(-^ production all of 

the final state is reconstructed. Also the value of this mass cut must be slightly higher 
than the equivalent for W'j (Eq. ( |5.6D ), to reflect the ~ 10 GeV higher mass of the Z 
boson. The lepton cuts are identical, without of course any requirement on the missing 
transverse energy. These cuts are motivated by an early CMS study [3]. Our results 
for ^/s = 7 TeV, shown in Table indicate that the Mg£^ cut is reasonably effective at 
removing the contribution to the cross section from photons in the Z decay. In the presence 
of the full lepton cuts, given in Eq. (|6.5| ), including photon radiation in the decay increases 
the cross section by about 15% at NLO. This is a larger difference than for VF7 production, 
which is to be expected since radiation may occur from both decay products. 

We now turn to the issue of the dependence of the cross section on the centre-of-mass 
energy y/s and the estimation of the theoretical uncertainty from scale variation. As is 
the case for the W'j cross sections of the previous section, we find that varying all scales 
by a factor of two about the central value of Mz results in a very small scale dependence. 
This is a result of the same accidental cancellation between the scaling behaviours of 
the component partonic cross sections. The scale uncertainties are therefore obtained by 
keeping Mp = Mz (since the fragmentation contribution is itself very small) and using 
{fiR = Mz/2,fj.p = 2Mz} and {fiR = 2Mz,fJ-F = Mz/2} for the upper and lower extrema 
respectively. Our results are shown in Table |5|. Finally, we can once again compare our 
NLO prediction for the total Zj cross section to a measurement already made at the LHC. 
From Table |^ we see that our NLO prediction for the cross section at 7 TeV and using the 
cuts of Eq. ( |6.3D is, 

a^^^iZj) X BR{Z ^ re+) = 9.83t[|-j5 pb . (6.6) 
The corresponding result reported by the CMS collaboration is [3], 

a'^^^{Z-f) X BR{Z = 9.3 ± 1.0 (stat) ± 0.6 (syst) ± 1.0 (lumi) pb (6.7) 
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Decay 


Cuts 


cj^O(e+e-7) 




No FSR 


Basic 7 


1.67(0) 


2.33(0) 




M^-y cut 


1.67(0) 


2.29(0) 




Lepton cuts 


0.82(0) 


1.17(0) 


Full 


Basic 7 


7.84 


9.83 




Mu^ cut 


2.08(0) 


2.81 




Lepton cuts 


0.99(0) 


1.39(0) 



Table 4: Z{-^ e^e")^ cross sections in picobarns at y^s = 7 TeV for the various scenarios described 
in detail in the text. Results in the upper half ("No FSR") correspond to neglecting diagrams 
containing photon radiation in the Z decay, while the cross sectio ns in the lower half ("FSR") 
include this effect. The cuts on the final state are specified in Eqs. (|6.3| )-(3.5). Statistical errors, 
unless otherwise indicated, are ±1 in the final digit. 



^/i [TeV] 


cj^O(e+e-7) [pb] 


a^^O(e+e-7) [pb] 


7 
8 
9 
10 
11 
12 
13 
14 


7.84(1) 
9.23(1) 
10.65(2) 
12.10(2) 
13.56(2) 
15.01(3) 
16.50(3) 
17.97(3) 


9.83(1)+^:?^: 

11.48(l)+^;^5 
13.10(l)+^.6% 

14.72(l)+^:^5 
16.38(2)+3-6% 

18.00(2)+3-5% 
19.61(2)+^-6% 
21.20(2)+3.7% 



Tab le 5 : Cross sections for Z(^- e+e^)7 production as a function of energy, using only the cuts of 
Eq. (|6.3[ ). The cross sections are calculated including the effects of photon radiation in the Z decay 
and the central values are obtained using /i^; — fip — Mp = Mz- The uncertainty is derived from 
the scale dependence, as described in the text. 



which is already in good agreement within errors. 

We conclude with an investigation of the importance of the gluon-gluon contribution in 
phenomenological studies. We shall use the full set of cuts given in Eq. (^^) as indicative of 
the appropriate experimental acceptance at the LHC. In that case we note that the relative 
effect of adding these diagrams is small since the cross section is dominated by regions of 
low that are enhanced for the qq process but not for the loop-induced gg diagrams. 
However, as one moves to moderate values of one would expect the relative size of the 
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Figure 9: The percentage of the Z{-^ e+e~)7 cross section using the cuts of Eq. < p.q ) contributed 
by the gluon-gluon initiated diagrams, as a function of the minimum photon px allowed in the 
events. The upper (red) curve is for ^/s = 14 TeV while the lower (blue) curve corresponds to 
y^=7 TcV. 



gluon-gluon contribution to grow. This is exactly the behaviour that we observe in Fig. |9|, 
with the gg fraction falling again at higher pj, due to the behaviour of the parton fluxes. 
We also see that, as expected, the gluon-gluon contribution is more important at 14 TeV, 
although it is still at most 3.5% of the total NLO cross section. 



7. WW production 

7.1 Description of the calculation 

The production of a pair of W bosons is an important channel, in part because of its role 
as a background to Higgs boson searches in which the Higgs decays into W pairs. The 
total cross section for the process, 

q + q^W+W- , (7.1) 

was first calculated in the Born approximation in ref. [40], with strong corrections to it 
given in refs. [41-43]. These processes are included in MCFM at NLO using the one- loop 
amplitudes presented in ref. [34]. Phenomenological NLO results for the Tevatron and the 
LHC operating at = 14 TeV have been presented in refs. [1, 44]. NLO results are 
also available for the processes W~^W~+ jet [45, 46], W"^W~^ + 2 jets [47] and W~^W~ + 
2 jets [48]. 

The contribution for the process, 

g + g^W+W-, (7.2) 
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[TeV] 


a^O{W+W-) [pb] 




7 
8 
9 
10 
11 
12 
13 
14 


29.51(1) 
35.56(1) 
41.75(2) 
48.07(2) 
54.53(2) 
61.10(3) 
67.74(3) 
74.48(4) 


47.04(2)+^JJ 

57.25(2)+^:^ 
67.82(3)+|8% 

78.70(3)+3;6j 

89.80(4)+3-3% 

101.14(5)+^-i% 

112.64(5)+3-0% 

124.31(6)+^-8% 



Table 6: Total cross sections for WW production as a function of energy. Renormalisation and 
factorisation scales are set to Mw- Upper and lower limits are obtained by varying the scales by a 
factor of two in each direction. Vector bosons are kept on-shell, with no branching ratios applied 

was first calculated in refs. [18, 19]. A more recent analysis of these contributions is given 
in ref. [21] where off-shell effects of the vector bosons and their subsequent decays are 
taken into account. Finally, the most complete analysis of these contributions to date is 
given in ref. [22] where the effect of massive quarks circulating in the loop is included. The 
authors find that the effect of including the third (t, b) isodoublet increases the gluon-gluon 
contribution by at most a factor of 12% at the 14 TeV LHC. 

The results presented in this section are obtained using our current implementation in 
MCFM which is as follows. Strong corrections to the qq initiated process are fully included, 
with additional contributions from singly resonant diagrams as described in ref. [1]. Since 
the contribution from the {t, b) isodoublet to the gg initiated process is small - certainly 
much smaller than the residual uncertainty resulting from the 0{a1) nature of the contri- 
bution - the gg process is included for two massless generations only. Our results for the 
gg process are in complete agreement with the equivalent two generation results presented 
in ref. [22]. As can be seen from Table 2 therein, the final cross section summed over qq 
and gg channels is smaller than the three generation result by 0.5%. 

The inclusion of the gg contribution with massless quarks in the loop is straightforward. 
The amplitudes can be obtained by simply recycling compact analytic expressions for 
certain contributions to the process e^e~ — )• 4 partons presented in ref. [27]. The precise 
relations are given in Appendix 

7.2 Results 

We begin our discussion of WW production by presenting the cross section as a function 
of ^/s in Table ^. The values are obtained by evaluating the cross section with a central 
scale choice of fiR = hf = M\y. Scale dependence is illustrated by presenting percentage 
deviations from the central value as the scales are changed simultaneously by a factor of 
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two in each direction. The W bosons are kept exactly on-shell and no decays are included 
for the cross sections presented in this table. We note that as for the other diboson cross 
sections the NLO corrections are typically large, enhancing the LO prediction by about 
a factor of 1.6. From the table, our NLO prediction for the total WW cross section at 
^ = 7 TeV is, 

^TVLO^ 47 0+2.0 (7.3) 

Although the general-purpose detectors at the LHC have collected only a handful of such 
events, both ATLAS [4] and CMS [5] have already reported first measurements of this cross 
section. They find, 

a^^^^^{WW) = 4ll^° (stat) ± 5 (syst) ± 1 (lumi) pb , (7.4) 
a'^^^^iWW) = 41.1 ± 15.3 (stat) ± 5.8 (syst) ± 4.5 (lumi) pb , (7.5) 

both of which are clearly compatible with the SM prediction. 

A measurement of the WW cross section at the LHC typically involves a jet veto to 
reduce the abundant top background [4, 5]. Since a jet- veto can change the relative size 
of the NLO corrections we will study the dependence of the NLO cross section on the 
transverse momentum scale used to veto jets, p'!^*'°. For our purposes here we define the 
jet veto as a veto on all jets with pT > p^^° that satisfy the rapidity requirement \r]j\ < 5. 
It is useful to consider the action of the jet veto under two sets of cuts, 

Basic WW : > 20 GeV, \r]^\ < 2.5, E^''' > 20 GeV , (7.6) 
Higgs : Basic WW + mu < 50 GeV, A(/)« < 60°, 

These cuts are typical of those used at the LHC to measure the total WW cross section, 
(with the additional application of a jet-veto) and those used to search for a Higgs boson [5]. 
The precise nature of the Higgs search cuts are dependent on the putative mass of the Higgs 
boson so here we have selected a set used for rrin = 160 GeV, when the decay to WW is 
largest. 

The ratio of the NLO to LO cross sections, as a function of p^^*° and for the two sets 
of cuts above, is shown in the upper panels of Fig. |l^. Since the gg initiated contribution 
does not contain any final state partons it is unaffected by the jet-veto at this order. As 
a result the relative importance of this contribution increases when a jet-veto is applied. 
We illustrate this by presenting the ratio a{gg)/a^^^ in the lower panels of Fig. |l^. As 
expected, the application of a jet-veto can reduce the iC-factor considerably. For instance, 
applying a jet veto at p1^^" = 20 GeV reduces the inclusive iC-factor by around 40%. From 
Fig. ^ we also observe that the Higgs cuts increase both the impact of NLO corrections and 
the gluon initiated contributions. The importance of the gluon initiated terms for Higgs 
searches has been observed in previous studies [26]. Indeed these studies have shown that, 
at ^/s = 14 TeV and with stricter cuts than those of Eq. (|7.7| ), the gg contributions can be 
as large as 30% of the NLO cross section [26] . At = 7 TeV and with cuts appropriate 
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Figure 10: The ratio of NLO to LO (upper) and the percentage of the NLO cross section from 
the gg initial state (lower) for WW — > e+/i~i^e£u production, as a function of the jet- veto pl^*". 
Results are shown using the basic cuts of Eq. ( |7.6| ) (upper, blue curves) and the Higgs search cuts 
of Eq. (7J) (lower, red curves). The NLO to LO ratio and gluon percentage with no veto applied 
are shown as dashed lines on the plot. The dashed lines are thus the asymptotic values of the solid 
curves. 



y/s [TeV] and cuts 


a^O(^e+fi-u,17^) [fb] 


a^^«(e+/.-z.eI7^) [fb] 


ii'-factor 


% gg 


7 (Basic) 


144 


249 


1.73 


3.05 


7 (Higgs) 


7.14 


15.19 


2.13 


6.85 


14 (Basic) 


296 


566 


1.91 


4.73 


14 (Higgs) 


13.7 


34.7 


2.53 


10.09 



Table 7: W+{- 



)W (— >■ (1 v^) cross sections in femtobarns at LO and NLO, the resulting 



if-factor and the percentage of the NLO cross section originating from gluon initiated contributions. 
Results are shown for the Basic (Eq. (7.6)) and Higgs (Eq. (7.7)) cuts. 



for this center of mass energy we find that the gg contribution is around 12% of the total 
NLO cross section with a jet veto of 20 GeV, as shown in the lower panel of Fig. 10. The 
values of the asymptotic limits of the jet veto curves shown in Fig |l^, corresponding to the 
ET-factor and gg percentage with no veto applied, are collected in Table |^. For completeness 
we also include the corresponding predictions for the cross sections at LO and NLO. 



- 21 - 



8. W^Z production 

8.1 Description of the calculation 

The production of a WZ pair proceeds at LO through the process, 

q + q'^W'^Z. (8.1) 

This process was first calculated to NLO in refs. [49, 50]. The inclusion of subsequent W 
and Z decays was added in ref. [43], partially including the effect of spin correlations. The 
full effect of spin correlations at NLO was later examined in refs. [1, 44], using the virtual 
amplitudes of ref. [34]. The QCD corrections to the process in which an additional jet is 
radiated are also now known [51]. 

The results presented in this section are obtained using the same implementation in 
MCFM as described in ref. [1]. In particular we include contributions from singly resonant 
diagrams that can be significant when one of the bosons is off-shell. The program includes 
both the contribution of a Z and a virtual photon, when considering the decay to charged 
leptons. We note that charge conservation precludes any contribution from gluon-gluon 
diagrams of the type previously discussed for WW production. 

8.2 Results 

The production of WZ pairs provides a valuable test of the triple gauge boson couplings 
(for a recent example, see for instance ref. [52]) and is a source of SM background events, 
for example in SUSY trilepton searches [53, 54]. It is also a background for SM Higgs 
searches in the case of leptonic decays, when one of the leptons is missed. In order to 
normalize the WZ background to such searches, in Table |8| we show results for the total 
cross section for WZ production at the LHC, as a function of the centre of mass energy. 
Both renormalisation and factorisation scales are set to the mean vector boson mass, {Mw-\- 
Mz)/2. Since the LHC is a proton-proton machine, the W~^Z and W~Z cross sections 
are not equal, with the ratio a^^'^ {W~ Z)/a^^^ {W~^ Z) varying between 0.56 (for ^/s = 
7 TeV) and 0.65 (at 14 TeV). 

The study of boosted objects at the LHC has potential as an additional handle on 
searches for new physics, such as a Higgs boson [55] or supersymmetric particles [56]. A 
possible first step for such searches would be to validate the method by performing a similar 
analysis for known Standard Model particles. In this regard, WZ production would be a 
natural proxy for associated Higgs production, WH, where the decay of the Z boson to a 
bottom quark pair is a stand-in for the decay of a light Higgs boson. 

To this end, in Fig. ^ we show the cross section for WZ production as a function of 
the minimum Z boson transverse momentum, at ^/s = 7 and ^/s = 14 TeV. We observe 
that the number of boosted Z bosons is sensitive both to the order in perturbation theory 
and the operating energy of the LHC. At both centre of mass energies the NLO prediction 
increases the number of high-p-p Z bosons, although the effect is larger at 14 TeV. To 
emulate a typical boosted Higgs search, with a px cut at 200 GeV, one thus retains about 
5% of the total NLO cross section - similar to the fraction for a putative Higgs signal. 
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[TeV] 


a^O(^W+Z) [pb] 


aNLO^^r+z) [pb] 


a^o^W-Z) [pb] 


^7VLO(p^-^) [pb] 


7 
8 
9 
10 
11 
12 
13 
14 


6.93(0) 
8.29(1) 
9.69(1) 
11.13(1) 
12.56(1) 
14.02(1) 
15.51(2) 
16.98(2) 


11.88(1)^^:^5 

14.48(l)+^.2% 

17-18(1)1^^5 
19.93(l)+^;85 

22.75(2)1^:^5 

25.63(2)1^:35 

28.55(2)13:2% 
31.50(3)1^:^5 


3.77(0) 
4.65(0) 
5.57(0) 
6.53(0) 
7.51(0) 
8.51(1) 
9.53(1) 
10.57(1) 


6.69(0)1^:^5 
8.40(0)1^:^5 

10.21(0)1^:^5 

12.11(1)1-5 
14.07(1)1-^5 
16.10(l)ltS 

i8.i9(i)l^:g 

20.32(1)1^:95 



Table 8: Total cross sections for WZ production as a function of energy. Renormalisation and 
factorisation scales are set equal to the average mass of the W and Z i.e. fJ.R = ^j.f = {Mw + Mz)/2. 
Upper and lower percentage deviations are obtained by varying the scales around the central scale 
by a factor of two. The vector boson are kept on-shell, with no decays included. 
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Figure 11: The fraction of the total WZ cross section surviving a cut on the Z boson transverse 
momentum, pt{Z) > p??'", at = 7 TeV (left panel) and 14 TeV (right panel). The NLO 

prediction is shown as a solid red curve and the LO one is dashed blue. 



9. ZZ production 

9.1 Description of the calculation 

Although the production of Z pairs is much smaller than the other diboson cross sections 
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considered above, it still plays an important role as principal background to searches for a 
Higgs boson around the Z pair threshold. The NLO corrections to the process, 



ZZ 



(9.1) 



were first calculated in refs. [57, 58], while the inclusion of spin correlations in the decays 
and phenomenology for the Tevatron and 14 TeV LHC was presented in refs. [1, 44]. 
Contributions from a gluon-gluon initial state, 



were first considered in refs. [18, 20]. The inclusion of leptonic decays of the Z bosons was 
examined in refs. [23, 24] and later investigated in the context of Higgs boson searches [25, 
26]. Furthermore, NLO results are also available for the closely-related ZZ+jet process [59]. 

The results presented in this section are obtained using our current implementation in 
MCFM which is as follows. Strong corrections to the qq initiated process include singly- 
resonant contributions - a slight extension of the results presented in ref. [1] - and the 
gg process is included for five massless fiavours. The contribution from massive top quark 
loops is suppressed by 1/m^ and is therefore neglected. This approximation results in 
gluon-gluon contributions that are 1% lower than those reported in refs. [25, 26], where 
the effects of massive top and bottom loops are included.^ Finally, we observe that all our 
amplitudes also contain contributions from virtual photons. 

The basic amplitudes entering the calculation of the gg contribution are simply related 
to those already discussed for the gg — ?■ WW process and are detailed in Appendix |^. 

9.2 Results 

We first present results for the dependence of the total cross section for ZZ production 
as a function of ^/s. As was the case for similar studies in previous sections we keep the 
Z bosons on-shell and do not include any decays. We choose a central scale choice of 
I^R = fJ-F = and vary this central scale by a factor of two in each direction to obtain 
an estimate of the theoretical uncertainty. Our results are shown in Table 

The decay of a Higgs boson to two Z's, which subsequently decay to leptons, is a 
promising search channel at the LHC. This is due to the fact that the Higgs will decay 
to Z's (with a moderate branching ratio) over a large range of Higgs masses that are not 
presently excluded. In addition, the four lepton signature associated with ZZ decay is 
experimentally clean. With Higgs searches in mind we apply the following cuts. 



In this definition of the cuts, ii and £2 represent the two hardest leptons and ^3 and 
^4 represent the two sub- leading leptons. The relevant distribution for the Higgs search is 
the invariant mass of the four- lepton system {m^i), for which we present our predictions in 
Fig. We show NLO predictions for both ^/s = 7 TeV and ^/s = 14 TeV, as well as the 



g + g^ZZ 



(9.2) 



p'r^'"^ > 20 GeV, > 5 GeV, \r]e\ < 2.5, mu,me'P > 5 GeV . 



(9.3) 



^We note that, when restricting our calculation to four massless flavours, our results are in complete 
agreement with the equivalent cross section quoted in ref. [25]. 



[TeV] 


a^O(ZZ) [pb] 


aNLO(^ZZ) [pb] 


7 
8 
9 
10 
11 
12 
13 
14 


4.17(0) 
5.06(0) 
5.98(0) 
6.93(0) 
7.90(0) 
8.89(1) 
9.89(1) 
10.92(1) 


6.46(0)+tS 

7.92(0)1S 
9.46(0)1-5 

11.03(0)1-5 

12.65(1)1S 

14.31(l)lf,5 
15.99(1)13.^% 

17.72(1)13;^% 



Table 9: Total cross sections for ZZ production as a function of energy. The renormalisation scale 
and factorisation scales are = = Mz- Vector bosons are produced exactly on-shell and no 
decays are included. 
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Figure 12: The invariant mass of the four lepton system in Z/j*{-^ e+e^)Z/7*(— > pro- 
duction at -y/s = 7 and y/s = 14 TeV, with the cuts of Eq. (|9.3D. In the upper panel we show both 
the total NLO prediction (upper curves) and the contribution from the gg initial state only (lower 
curves). In the lower panel we plot the fraction of the NLO prediction resulting from the gg initial 
state. 



contribution from the gluon-gluon diagrams alone. 

From the figure we observe that, although the gluon initiated pieces are fairly important 



at the level of the total cross section, their effect in the region m^i < 2Mz is rather smaller 
(at the few percent level). As this threshold is crossed the percentage effect increases to 
around 7% (7 TeV) or 10% (14 TeV). Our results at 14 TeV agree with the findings of a 
previous study in a similar kinematic range [26]. It is clear that the gg initiated piece is 
most important as a background to Higgs bosons searches in the region niH > 2Mz- 

10. Conclusions 

In this paper we have provided NLO predictions for all diboson processes at the LHC, 
both at the current operating energy of ^/s = 7 TeV and at higher energies appropriate 
for future running. The calculations are contained in the parton level code MCFM, which 
includes the implementation of pp — t- 77 for the first time. In addition, where appropriate 
we have revisited the treatment of many of the vector boson pair processes in order to 
ensure the relevance of the predictions for the LHC. 

In order to enable simpler comparisons with experimental results we have implemented 
experimental photon isolation cuts into MCFM. This requires the inclusion of fragmenta- 
tion contributions [6, 8], in which a QCD parton fragments into a photon plus hadronic 
energy. These fragmentation contributions require the introduction of fragmentation func- 
tions that contain both non-perturbative and perturbative information. Including this 
isolation condition extends the previous treatment of photons in MCFM, for which the 
smooth cone isolation of Frixione [9] had been used. Although this latter method is simple 
to implement theoretically it is not well-suited to experimental studies. 

At the LHC, contributions to diboson production which proceed through a gluon ini- 
tiated quark loop can have a significant effect on cross-sections. Although formally in 
perturbation theory they enter at NNLO the large flux of gluons at LHC center of mass 
energies can overcome the formal 0{q^) suppression. Consequently we have included the 
gluon initiated processes gg — >■ {77, Z^, ZZ, WW} whose contributions have been stud- 
ied in the past [11-26]. We have also included higher order corrections to the gg 77 
process [13, 14], which are formally at the level of N^LO in perturbation theory. These 
corrections are indeed of phenomenological relevance at the LHC, since at "NLO" the gluon 
contribution is around 20% of the total cross section at ^/s = 7 TeV. 

We have presented detailed results for the diphoton process, pp — > 77 at the LHC, 
which is an extremely important channel for light Higgs boson searches. We have presented 
theoretical predictions for the rrijj distribution using experimental cuts and isolation. We 
have illustrated how these cuts rediicc the nominal cross section for a range of \/s appro- 
priate to the LHC and also investigated the sensitivity of our predictions to the amount of 
hadronic energy in a fixed size isolation cone. We have also shown that experimental Higgs 
search cuts, which usually require staggered photon transverse momenta, produce large 
if-factors at NLO due to the limited kinematic configurations probed at leading order. 

We also presented results for Wj and Z7 production. As a result of the new isolation 
procedures we were able to compare our NLO prediction for the cross sections with the 
recently measured values from CMS. We investigated the effects of various lepton cuts on 
cross-sections and predictions for distributions, particularly those cuts designed to suppress 
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the contribution of photon radiation in the vector boson decay. This is an important 
consideration in the search for anomalous coupUngs between vector bosons and photons. 
Although the gluon-gluon contribution to the cross section is a few percent of the total 
at NLO, we found that as the minimum photon pT is increased the gluon initiated terms 
become relatively more important. 

We studied the effects of including gluon initiated processes on WW and ZZ produc- 
tion. These pieces have been calculated using various methods in the past [21, 22, 25, 26], 
where it has been shown that the effects of the massive top quark are small. For this rea- 
son we have ignored the effect of a third generation (for WW) or top quark loop (for ZZ) 
and instead only include loops of massless quarks using the analytic formulae described in 
Appendix ^ For ZZ production we also include the effects of singly-resonant diagrams 
that had previously been neglected in ref. [1]. For the case of WW production we paid 
particular attention to the effects of a jet-veto on the NLO cross section, since at the LHC 
a jet-veto is necessary in order to reduce the abundant top background. We found that 
applying a jet- veto and Higgs search cuts increases the overall percentage of the cross sec- 
tion associated with the gluon initiated process. For ZZ production we found that the 
contribution of the gg process is only important in searches for a Higgs boson with a mass 
greater than 2Mz- The fact that these gg corrections can be large in some circumstances, 
suggests that the two loop corrections are worth calculating, to get a better idea of the 
associated theoretical error. 

We also presented results for WZ production at the LHC. As an example, we investi- 
gated the fraction of events that survive a minimum cut on the Z transverse momentum. 
This quantity is important in regards to boosted searches for Higgs bosons and supersym- 
metry. 
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A. Input parameters for phenomenological results 

The electroweak parameters that we regard as inputs are. 



Using the values of M\y, Mz and as above then determines Oe.m.l-^z) and sin^ as 
outputs, where 6^, is the Weinberg angle. We find. 



Mw = 80.398 GeV , Mz = 9L1876 GeV 
Vw = 2.1054 GeV , Tz = 2.4952 GeV , 
Gf = 1.16639 X lO-^GeV-2 . 



(A.l) 
(A.2) 
(A.3) 



sin^ (9^ = 1 - M'^/MI = 0.222646 



(A.4) 
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= = 132:338 ■ (^-'^ 

This value of ae.m. may not correspond to the value of Oe.m. used to fit the fragmentation 
functions in ref. [7] . The value of ae.m. used in their fit is hard to extract from ref. [7] . After 
isolation we believe that any potential mismatch will be of minor numerical significance. 

B. Helicity amplitudes for gluon-gluon processes 

In this appendix we present results for three of the gluon initiated processes considered 
in the text, namely gg — )• Z7, gg — )• WW and gg — )• ZZ. We first describe some general 
notation and then consider each of these processes in turn. 

B.l Notation 

In order to specify the amplitudes we first introduce some notation. The QED and QCD 
couplings are denoted by e and gs respectively and Q'^ is the charge of quark q in units of 
e. The ratio of vector boson V (= W, Z) and photon propagators is given by, 

s - Mp + I Tv My 

where My and Ty are the mass and width of the boson V. Fermions interact with the Z 
boson through the following left- and right-handed couplings, 

-l + 2sin2 0^ 2sin2^^ 

'"L = — T^Ti , Vr - 



sin26',„ ' " sin 261, 



w 



, _ ±l-2Q'?sin^ g _ IQisin^e^ 

sin 20^ ' sin 20^ ' ^^'^^ 

The subscripts L and R refer to whether the particle to which the Z couples is left- or 
right-handed and the two signs in correspond to up (-I-) and down (— ) type quarks 
We express the amplitudes in terms of spinor products defined as, 

(ij) = u^{pi)u+{pj), [ij] = u+{pi)u^{pj), (ij) [ji] = 2pi ■ pj . (B.3) 

B.2 Amplitudes for gg — )• Z7 

In this section we present results for the amplitudes relevant for the process, 

^ g{pi) + giP2) + 7(^3) + i{p4) + i{P5) ■ (B.4) 

These amplitudes can be extracted from the fermion loop amplitudes ^g'"^(lq, 2^, 3^,4^, 5j, 6^) 
in ref. [27], by taking the limit in which the quark and antiquark are collinear, and renaming 
the momenta. The boxes and triangles which potentially could contribute to this ampli- 
tude are shown in Fig. |l^. Considering the triangle diagrams first, we see that diagrams 
like Fig. |l3|(c) can never give a contribution because of colour. Fig. |l3|(b) with the vector 
coupling of the Z to the fermion loop vanishes because of Furry's theorem, and with the 
axial coupling to the fermion loop vanishes because of Bose statistics (Landau- Yang the- 
orem). We therefore only have to consider box diagrams like Fig. |l^(a). In this case the 
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(a) 




(b) 



(c) 



Figure 13: Examples of diagrams that could potentially contribute to gg — > ^U. 



diagrams with an axial coupling vanish, so we only have to consider box diagrams with a 
vector coupling. 

The result for the fully dressed amplitude is, 

^2 



nf 



1=1 



167r2 

-Q' + lvlR{vi + v'n)rz{s,G) 



(B.5) 



where ai,a2 are the colour labels of the two gluons and there are Uf flavours of massless 
quarks circulating in the loops. 

The amplitude A^\l^ ,2^ ,3'^ ,4:J is entirely rational and given by the following 
expression, 

^^(i+,2+,3+,47,5+) = 2^ , r." - ^ , +1^2 



(12)' (13) (4 5) 2 (12)' [5 4] 



(B.6) 



The amplitude A^^\l'^ , 2+, 3^ , 4^ , 5^) contains dependence on the box and triangle func- 
tions Lo(r), Li(r) and Ls_i(ri,r2) that will be defined below. The result is. 



+ 



(12)^(45) 
'^ (2 3) (14) (2 4) [21] ^ /-S123 

Z — -_ — — Lq 



[31] (12)^ (4 5) 
(13) (2 4) [21] [51] - 
[31] (12)' (4 5) [5 4]. 



+ 



-Sl3 

1 o 2 



•S123 — •S123 

(2 4)' [2 1]' 
[3 1]'(12)'(4 5) 



Li 



-^123 

(B.7) 
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We note that the 1 -f-)- 2 symmetry in this equation is to be apphed to the terms inside 
square brackets only. 

The final helicity amplitude A^^'^ ( 1+ , 2^ , 3:^ , 4^ , 5^ ) can be obtained by exchange from 



Eq. dBj), 

4'')(l+,2-,37,47,5+) 



J [3lf[25]^ + ]35]W 
I [3 2]^ [5 4] 

' J21] [3 5] [2 5] (2 3) ^ 
(13) [3 2]^ [5 4] ° 
[31] [2 5] (2 3) (43) 
(13) [3 2]2[5 4] (4 5) 



Ls_ 



-Sl3 



-Sl2 



+ 



-gl23 
-Sl3 

2^3 



-Sl23 — Sl23 

[2 5]^ (2 3)^ 
(13)^ [3 2]^ [5 4] 



Li 



-^123 
-Sl3 

(B. 



where again the 2 -f-)- 3 is to be applied to the terms inside square brackets only. 

The latter two amplitudes are defined in terms of the following functions that arise 
from box integrals with one non-lightlike external line, 



Lo(r) 
Li(r) 



ln(r) 
1 — r 
Lo(r-) + 1 
1 — r 



vr 



Ls_i(ri,r2) = Li2(l - ri) + Li2(l - + Inri lnr2 — 

b 



and the dilogarithm is defined by, 

Li2(x) 







y 



(B.9) 
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The remaining amplitudes can be obtained from these ones by simple symmetry operations. 
B.3 Amplitudes for gg — )• WW 

In this section we present results for the amplitudes relevant for the process, 

^ g{pi) + g{P2) + Mps) + Apa) + ^'(Ps) + ^e'iPe) ■ (B.ll) 



There are six contributing Feynman diagrams depicted in Fig. 14. These diagrams rep- 
resent exactly the same set that appears in the calculation of certain contributions to the 
process e~^e~ — )• 4 partons. We therefore simply reinterpret the compact expressions for 
such amplitudes presented in ref. [27], modifying the overall factor appropriately. Specifi- 
cally, we find that the contribution from a single generation of massless quarks in the loop 
is given by. 



-i-loop/,hi o- 4+ 5- a+ 



j;aia2 / 9w9s 

167r2 



rw{s3A)rw (^56) ^6:4 ( 3+, 4r , l"'^! , l'''' ; 6r , 5, 



. (B.12) 
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Figure 14: Diagrams that contribute to gg — >■ WW. 



The helicities and colour labels of the two gluons are hi, /12 and ai, 02 respectively and the 
amplitude Aq.^ is defined in Sections 2 and 11 of ref. [27]. The operation [. . .^'^ exchanges 
the spinor products (a 6) and [ab] in the amplitude. The labelling on the right hand side 
of this equation is as written in ref. [27]. For our purposes we make the identification on 
the left hand side, ((7 — t- g — )• ^, e — ?■ i' and e — ?■ I^^/). 

B.4 Amplitudes for gg — > ZZ 

In this section we present results for the amplitudes relevant for the process, 

^ g{pi) + g{P2) + iiPs) + e{pi) + e'{p5) + I{P6) ■ (B.13) 

The extension of the procedure outlined above for gg — )■ WW is clear. One must now 
simply sum over all four possible helicity combinations for the leptonic decays. The result 
for the fully dressed amplitude is, for a particular choice of lepton helicities. 



^1— loop 

5' 



iftl r,/l2 o- A+ C + 



4 2 

aia2 ^ 9s ^ ^ 



27r2 



i=l 



Q' + i^vIr{vI + v'r)Vz{sm) 



-Q' + i^viR{'"l + v\)Vz{s,^) 



0+ A- 1-/11 9-/12. c- c; + 



(B.14) 



The labelling on the right hand side of this equation is as written in ref. [27]. For our 
purposes we make the identification on the left hand side, {q ^ £, q ^ i, e ^ £' and 
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